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Abstract. We present a study of the rnr scattering amplitude in the a and p channels at finite tempera- 
ture and nuclear density within a chiral unitary framework. Meson resonances are dynamically generated 
in our approach, which allows us to analyze the behavior of their associated scattering poles when the 
system is driven towards chiral symmetry restoration. Medium effects are incorporated in three ways: 
(a) by thermal corrections of the unitarized scattering amplitudes, (b) by finite nuclear density effects 
associated to a renormalization of the pion decay constant, and complementarily (c) by extending our 
calculation of the scalar-isoscalar channel to account for finite nuclear density and temperature effects in 
a microscopic many-body implementation of pion dynamics. Our results are discussed in connection with 
several phenomenological aspects relevant for nuclear matter and Heavy-Ion Collision experiments, such 
as p mass scaling vs broadening from dilepton spectra and chiral restoration signals in the a channel. We 
also elaborate on the molecular nature of tttv resonances. 



PACS. Il.lO.Wx Finite temperature field theory 
- 25.75.-q Relativistic heavy-ion collisions 



12.39.Fe chiral lagrangians - 21.65.+f Nuclear matter 



1 Introduction. 

The lightest meson resonances, the p{77Q) and the /o(600) 
or (T, play a crucial role in different phenomena pertaining 
to the hot and/or dense medium created both in Rela- 
tivistic Heavy Ion Collisions and in Nuclear Matter exper- 
iments. The modifications of the spectral function of the p 
resonance in medium are crucial to understand correctly 
the dilepton yield emerged from Heavy Ion Collisions [HIS] . 
The two main theoretical scenarios currently proposed in 
the literature can be classified into resonance mass shifting 
and broadening, according to the resulting dominant ef- 
fect. Mass shifting models are inspired by the Brown-Rho 
(BR) scaling hypothesis [3 , which predicted that vector 
meson masses should scale with the quark condensate and 
therefore the main spectral modification of the resonance 
would be dictated by chiral symmetry restoration. This 
scenario is supported also by the so called hidden local 
symmetry approach \4^. The broadening-dominated sce- 
nario is supported by different theoretical analyses [SUHllIl 
[8l[9] including recent Unitarized Chiral approaches [IOllTT| 
I12j . The most recent experimental dimuon data from the 
NA60 Collaboration 2 clearly favor a broadening situa- 
tion with negligible mass shift, whereas the earlier CERES 
results [1] were reasonably explained by both descriptions. 
It is worth mentioning also the results of the STAR Col- 
laboration at RHIC [13 , which has reported a sizable mass 
reduction by medium effects measured in p*^ — > tt+tt^ in- 
stead of dileptons. The modifications of the p properties 



Electronic address: gomez@fis.ucm.es 



have been also measured in cold nuclear matter experi- 
ments. The E325-KEK collaboration 14J has reported a 
measurable shift in the masses of vector mesons compati- 
ble with theoretical predictions based on Brown-Rho scal- 
ing |TS] and QCD sum rules [IS]. On the other hand, the 
JLab-CLAS experiment [17] has obtained results compat- 
ible with vanishing mass shift, as predicted by most of 
in-medium hadronic many-body analyses where broaden- 
ing is the dominant effect [S l lTllTO l lTS] . 

The possible modification of the /o(600)/(t in hot and 
dense matter is interesting because this is a state with the 
same quantum numbers as the vacuum and therefore it 
might be sensitive to chiral symmetry restoration. In this 
sense, an early proposal [T9] suggested that the a could in- 
duce a measurable threshold enhancement of the tttt cross 
section, which would be interpreted as a precursor of chiral 
symmetry restoration. The argument was that the mass of 
the a state should decrease by medium effects, since it is 
proportional to its vacuum expectation value in the chiral 
limit. Such a decrease would eventually shrink the avail- 
able two-pion phase space when the a mass reaches the 
two-pion threshold, producing a bump in the imaginary 
part of the scattering amplitude, due to the proximity of 
the pole to the real axis. It is important to remark that in 
this original argument, it is implicitly assumed that i) the 
cr is dominated by its qq component so that its expecta- 
tion value behaves like the quark condensate and ii) that 
the a is narrow enough so that its width (imaginary part 
of the pole) vanishes when its mass (real part) approaches 
the threshold. None of these assumptions seem to be sup- 
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ported by the physical (vacuum) /o(600) state quoted by 
the Particle Data Group [20 , which is a very broad state 
measured in tttt scattering. Although it is commonly ac- 
cepted that this state is a member of the scalar nonet, its 
qq nature has been criticized on the basis of lattice [H] and 
large- TVc [22] analyses. The physical state is likely to have 
important won-qq compound such as tetraquark, glueball 
or meson-meson, commonly referred to as a "molecular" 
state. Nevertheless, threshold enhancement is indeed ob- 
served in nuclear matter experiments, both in it A — s- tttt A' 
and in 7 A tttt A' gS] reactions. Although the 
size of the effect is still under debate, a clear signal is seen 
in the scalar channel for increasing nuclear density as com- 
pared with the vector channel and this is in fair agreement 
with most theoretical analyses at finite nuclear density [26l 
[?7l[^8. 29 , 30,.,31j . In contrast, finite temperature analyses 
show that this state remains broad even near the chiral 
phase transition despite the proximity of the pole to the 
two-pion threshold [12,30,32J, which in practice does not 
produce any sizable enhancement in the scattering ampli- 
tude or cross section. A striking possibility to be explored 
is that by increasing further the medium strength, the a 
could become a tttt bound state, as suggested earlier in 
[33] and confirmed recently in [12,.30,,32J. 

In the present work, we will investigate further about 
these issues, within the context of unitarized chiral ap- 
proaches. The main goal is to establish to what extent 
chiral symmetry dictates the in-medium properties of the 
light meson resonances. The most general framework to 
account for all the interactions compatible with chiral 
symmetry is the effective chiral lagrangian approach. The 
most prominent example is Chiral Perturbation Theory 
(ChPT) for the meson sector [34,35j, but it can be equally 
applied to the meson-baryon one [36]. Since these effec- 
tive theories are built basically as expansions in deriva- 
tives or energies, they cannot account for resonances, since 
the chiral expansion violates the unitarity bounds. This 
has been traditionally solved by introducing unitarization 
methods, giving rise to the so called chiral unitary ap- 
proaches, which have proved to be very successful in vac- 
uum to describe meson-meson and meson-baryon interac- 
tions and generate dynamically low lying resonances [39|, 
I40j . Furthermore, as commented above, the unitarization 
program has been extended to account for finite temper- 
ature and density effects [IO l fTT ] fT2 l l29 l [3T ] [37 l l38] . 

Here we will analyze some of our recent results for the p 
and <T mesons obtained within the unitarized chiral frame- 
work, paying special attention to their nature and their 
role in chiral symmetry restoration and studying some 
aspects not considered before like a new analysis of the 
combined effects of temperature and nuclear density for 
the a meson based on a Lippmann-Schwinger (or Bethe- 
Salpeter) equation approach accounting for many-body 
pion dynamics versus a simplified /7r-scaling scenario, a 
direct comparison with nuclear matter experiments of our 
results for the p mass linear density dependence and the 
interpretation of the results obtained for the behavior of 
the resonances near threshold in terms of a "molecular" 
classification of those states. 



The paper is organized as follows: in Section [2] we will 
present the formalism and results within the framework 
of the Inverse Amplitude Method (lAM) at finite tem- 
perature. In Section [3] we will introduce nuclear density 
effects only by rescaling properly the pion decay constant 
at T = within the lAM. In that section, we will pro- 
vide in particular an interpretation of our results in terms 
of "molecular" classification and a numerical comparison 
with experimental results from dilepton decays in reso- 
nance production in finite nuclei. Finally, in Section |4] we 
present a new calculation for the a channel which includes 
temperature and nuclear density many-body effects in or- 
der to compare our different approaches. 



2 Finite temperature resonances with the 
Inverse Amplitude Method. 

One of the simplest and more powerful unitarization meth- 
ods for chiral theories is the so called Inverse Amplitude 
Method (lAM) [39'. Its name comes from the simple ob- 
servation that unitarity implies that the inverse of a given 
partial wave amplitude t^'^ in vrvr tttt scattering should 
satisfy: 

S^S^l^ Imt"(s) = CTo(s)K^^(s)|2 



for s > 4TOj,where s is the center of mass energy squared 
and ao{s) = -y/l — Am^/s is the two-pion phase space. 

Consider now the ChPT expansion of partial wave am- 
plitudes: 

t"{s)=4-'{s)+ti-'{s)+0{p'') (2) 

Here, p denotes generically a meson momentum, mass 
or temperature {p is to be compared with the character- 
istic chiral scale ^ 1 GeV, whereas T is meant to be 
below Tc ~ 200 MeV) and tk is the 0{p'') contribution. 
Recall that, according to the standard ChPT power count- 
ing [341135] . t2 accounts for tree level diagrams from the 
lowest order lagrangian €2 - Up to that order, only the pion 
decay constant /^r and the pion mass TOtt enter the result. 
The order includes the one-loop diagrams from £2 plus 
the tree level £4 terms needed for renormalization. The 
£4 low-energy constants entering the pion scattering am- 
plitude, when it is expressed in terms of the physical ttItt, 

are denoted li — in the convention of [35j . 

The chiral expansion ([2]) satisfies only a perturbative 
version of the unitarity relation ([1]), namely: 

lmti'is)^ao{s)\ti'is)\' (3) 

and so on for higher orders, which eventually means that 
chiral expansions are not compatible with the bounds on 
partial waves implied by unitarity. In other words, they 
grow arbitrarily with energy. Unitarization methods allow 
to construct chiral amplitudes that are exactly unitary. 
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In particular, the lAM amplitudes are built by demand- 
ing i) exact unitarity and ii) that at low energies they 
match the ChPT series to a given order. These conditions 
lead to the lAM result, which is formally justified by the 
use of dispersion relations [3^. There is however a fur- 
ther, more technical requirement, which is that the lAM 
partial-wave aplitudes should vanish at the same values 
of the energy and with the same power as the perturba- 
tive amplitudes. These values are the so called Adler zeros 
and lie below threshold. Since a zero of the amplitude is 
a pole of its inverse, this affects the analytic structure of 
1/t. A detailed discussion can be found in [41] where the 
proper correction to the lAM is derived using dispersion 
relations and it is shown that these additional terms pro- 
duce a negligible effect in the physical region. However, 
as discussed in |12j, taking into account this correction 
is important when dealing with medium effects that can 
drive the poles to the real axis, as it is the case here, since 
otherwise there would be spurious poles both in the first 
and second Riemann sheets below threshold. 

The lAM can be extended at finite temperature by 
including the thermal corrections to the scattering ampli- 
tude, which have been calculated in [42] to one loop in 
ChPT. Since temperature enters only in the loops, t2 is 
T- independent. For t4{s;T) one gets a perturbative uni- 
tarity relation exactly like ([2]) but with the phase space 
replaced by: 

ctt(s) =ao(s)[l + 2ns(\/i/2)] (4) 

with nsix) = [exp(x/T) — 1]^^ the Bose-Einstein distri- 
bution function. 

The function (tt(s) is the thermal phase space, which 
is increased with respect to the T = one by the dif- 
ference [1 + nB{Ei)] [1 + nBiE2)] - nBiEi)nB{E2) = 1 + 
nB{Ei) + nB{E2), where Ei^2 are the energies of the two 
colliding pions, corresponding to the difference between 
enhancement due to the increase of two-pion outgoing 
states and absorption due to collisions of the incoming 
pions with the thermal bath ones. In the center of mass 
frame, where partial waves are defined, Ei — E2 = 
and the thermal phase space reduces to One can then 
use the same T = lAM requirements, replacing ctq (Jt 
and the partial waves by the finite-T ones, provided that 
only intermediate two-pion states are relevant in the ther- 
mal bath, as expected in a dilute gas regime at low and 
moderate temperatures. Finally, one arrives to the ther- 
mal lAM formula for a given partial wave: 

,IAM ^ hjsY 

t2{s) - U{s-T) + A{s-T) 

A{s;T) = U{S2;T) - " ^^)(^ " _ ^^(,,.7)] 

S- SA 

(5) 

where the A function is the Adler zero contribution dis- 
cussed above, sa denoting the Adler zero (T-dependent) 
expanded as sa = S2+S4 + . . . with 54 = — t4(s2; T)/t'2is2)- 



Performing the conventional extension of the ampli- 
tude to the second Riemann sheet, one finds poles in the 
I = J = and I — J = 1 channels which are identi- 
fied as the /o(600) and the p{770). We show in Figure [1] 
the results for the pole position Spoie = {Mp — iFp/i)'^ 
for different temperatures. The li values we have used are 
li = —0.3, I2 = 5.6, I3 = 3.4 and I4 = 4.3, which give 
for the mass and width of the p(770) at T = Mp ~ 756 
MeV and Ep ~ 151 MeV. For the /o(600)/cr at T = we 
find Mp ~ 441 MeV and Fp ~ 464 MeV. 

The general features we observe are that the thermal 
p pole shows a predominant and increasing broadening 
behavior, while for the a an important mass decrease 
takes place, presumably due to chiral restoration, while 
the width increases for low temperatures but decreases 
for temperatures of T ~ 100 MeV and beyond. In the rest 
of this section, we will discuss in more detail these dif- 
ferent behaviors in connection with the phenomenological 
issues commented in the introduction. 



2.1 The thermal p meson: broadening versus mass 
scaling in tttt scattering and dilepton probes. 

The p pole obtained in our lAM thermal approach under- 
goes a significant broadening at finite temperature. The 
main source of thermal broadening is the Bose-Einstein 
increase of phase space given in However, it is not 
the only one. In fact, using the Breit-Wigner {Fp ^ Mp) 
parametrization for the p exchange in tttt —> p —> tttt, one 
gets [11]: 

^ = [l + 2n.(M./2)]f^ (6) 

where g is the effective ptttt vertex. 

In Figure [5] (left) we show in detail the dependence of 
the mass, width and effective vertex and is clearly seen 
that the width increases roughly as the phase space up 
to T ~ 100 MeV and then further broadening arises due 
to the increasing effective pir-K vertex. The broadening we 
have obtained is also present in the electromagnetic pion 
form factor [43], which enters directly into the dilepton 
yield (for back to back dileptons in our case) arising from 
7r~''7r~ . Our result is therefore compatible with the 

broadening scenario observed clearly in NA60 [2] (see the 
discussion in the introduction). 

The mass of the pole barely changes with temperature 
in our approach. This seems to be in contradiction with 
the dropping-mass scenario obtained in scaling models, 
which predict that the mass should scale roughly with the 
quark condensate Mt/Mq ~ {qq)T/{<iq)o [5]- However, it 
is worth mentioning that more recent analyses based on 
the same scaling hypothesis [4ll44] suggest that the mass 
dropping might be really effective only very close to the 
transition temperature, which in practice would mean that 
those predictions are not strictly incompatible with the 
dilepton data. 

In order to clarify our results on this matter, we have 
plotted in Figure [U (right) the pole mass and the quark 
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Fig. 1. Temperature dependence of the /o(600) and p complex poles unitarized by the lAM. The points are obtained by varying 
the temperature in 20 MeV intervals. 
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Fig. 2. Thermal width and mass of the p as extracted from its lAM pole. In the left panel we also show the T-dependence of 
the effective pirn vertex. In the right panel we compare the mass with the quark condensate. The mass vanishes at Tb ~ 310 
MeV. 



condensate calculated from a virial expansion f5S] using 
the 0{p'^) pion scattering amplitudes, which gives a crit- 
ical temperature Tc ~ 230 MeV. The results are extrapo- 
lated up to the temperature To where the mass vanishes, 
although our ChPT-based approach is not meant to be 
vaUd there. We see that the mass drops rather abruptly, 
with a power, while, as said before, it remains almost 
constant for T below the chiral transition. We do not see 
a scaling pattern when compared to the condensate and, 
besides, Tq ~ 310 MeV lies far from the critical value 
where the condensate vanishes. On the other hand, in BR- 
like scaling models, the effective vertex decreases and, 
as commented before, there is no significant broadening. 
In conclusion, although we obtain a dropping mass qual- 
itatively compatible with BR-like models, our dominant 
broadening effect, the increase of the effective vertex and 
the departure from the condensate are in conflict (at finite 
temperature) with that scenario. 



2.2 The thermal /o(600)/(t meson: threshold behavior 
and qq nature. 

As mentioned in the introduction, the interest in the a 
pole concerns mainly its role as a precursor of the phase 
transition. The temperature dependence we obtain in Fig- 
ure [1] shows that mass decreasing is a prominent feature 
in this channel, thus signaling chiral restoration. In fact, 
this is the dominant effect here over phase space increase 
for T ~ 100 MeV and beyond. However, and this is cru- 
cial as far as observable effects are concerned, the pole 
remains broad even when its real part has reached the two- 
pion threshold. The reason behind this is that for a broad 
state {Fp Mp) the usual result T cx a{]VP)e{AP -Ami), 
which forbids the pole to be at threshold with a nonzero 
width, should be replaced by F ^ f a{s)9{s — 4m^)p(s) 
[12j where p{s) is the resonance spectral function, different 
in general from the narrow (5-function p{s) = 27rS{s~ M'^). 
Therefore, the width can remain sizable near threshold. In 
fact, a simple model of thermal width based on this ob- 
servation turns out to reproduce quite well the obtained 
lAM poles [12]. This behavior for the a pole implies no 
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threshold enhancement at finite temperature as a pre- 
em'sor of the transition, since the pole is still far from 
the real axis. This is clearly observed when we plot the 
squared modulus of the I = J = partial wave in Fig- 
ure [S] (left). Threshold enhancement for a typical nar- 
row resonance would mean for the amplitude Imi(s) ^ 
2771^6' (s — 4to^) / ^/s — 4mJ for s^4mj [H]. As we shah 
see below, the situation changes dramatically when finite 
density effects are included. 

The role of the a pole as a precursor of the transi- 
tion is strongly linked to its qq nature, as discussed in 
the introduction. In Figure [3] we compare the mass of the 
pole with the root of the quark condensate, similarly as 
we did with the p in the previous section. Recall that, 
in a 0(4) model, ~ (cr) ~ /,r and, on the other 



hand, 



J (qq) I from the Gell-Mann-Oakes Ren- 



ner (GOR) relation [46]. We also plot the results near 
the chiral limit, where the condensate vanishes at a lower 
temperature and explicit chiral symmetry breaking effects 
are minimized. In any case, we do not see a scaling pat- 
tern, as expected from our result of a broad thermal state. 
Near the chiral limit, the pole mass does not even go to 
threshold near the critical temperature, since the width is 
notably increased due to the more available phase space. 
From our analysis of thermal effects one can also con- 
clude that the non-gg component of the /o(600) must be 
of crucial importance. As discussed in the introduction, 
the same conclusion has been reached in vacuum studies 
pT1[22]. The novelty here is the use of thermal arguments. 
We will see in Section [3] that finite density effects driving 
the resonance poles to the real axis below threshold al- 
low also to extract interesting conclusions about their qq 
nature. 



3 Finite density chiral restoring effects at 
T = through scaling in the I AM. 

The simplest way to incorporate nuclear density effects at 
T = in the unitarized chiral approach is to encode them 
only in the variation of the pion decay constant, to linear 
order in density, as [37]: 



!l{p) ^ {qq){p) 

fm ~ {mm 



1 - 0.35- 



Po 



where p is the nuclear density, a^^N — 45 MeV is the pion- 
nucleon sigma term and po ~ 0.17 fm"'^ is the normal or 
saturation nuclear matter density. 

Therefore, by varying in our lAM approach, we 
mimic chiral restoring nuclear effects [12] . This approach 
clearly ignores standard many-body corrections, like the 
coupling of pions to particle-hole {ph) and Delta-hole (Ah) 
excitations, considered for instance in [51 fTUl[^51l5T] . This 
approximation is meant to be more adequate for the a 
than for the p since, as we have seen at finite tempera- 
ture, chiral restoration tends to dominate the a pole be- 
havior under medium effects. In fact, by changing only 



/tt, no medium broadening is produced, which is probably 
unrealistic for the p case, as emphasized in many-body 
works [5l[7lll0|. This must be borne in mind when inter- 
preting our results for the p in terms of BR-scaling. In any 
case, in SectionUjwe will consider a different unitarization 
scheme, which allows to introduce all the above mentioned 
nuclear many-body effects and compare with our simple 
"/tt scaling" considered here. 

The results obtained by varying in the lAM ampli- 
tudes are displayed in Figs. |4| and [5| for the pole trajecto- 
ries in the complex plane and mass scaling, respectively, 
in the 00- and 11-channels. In Fig. [6| (right) we also show 
the effect on the scattering amplitude in the a channel. As 

decreases the a pole becomes narrow enough so that 
chiral mass reduction brings it to the real axis, which pro- 
duces threshold enhancement in the amplitude although 
at densities well above po. Regarding the mass dropping, 
it takes place now along with the condensate (see Fig. [5|) . 
In Fig. [7|we show |iooP considering simultaneously tem- 
perature and finite density according to Eq. ((7|) for f-^ . We 
observe that, for a given value of /tt (small enough), the 
net effect of introducing temperature is to amplify thresh- 
old enhancement, as the amplitude is notably softened at 
higher energies, by thermal broadening, what makes the 
low energy region relatively more important. Still, one has 
to keep in mind that (a) this effect is driven by the prox- 
imity of the (T pole to the real axis when decreasing 
(b) thermal effects alone do not generate threshold en- 
hancement in the amplitude as discussed above, and (c) 
the TTTT phase space is not open below 2m7r even at finite 
temperature. In Sect. [4] we will compare these results for 
the 00-channel with an implementation of nuclear density 
effects in a dynamical many-body calculation. 

As for the p channel, the BR-like scaling pattern is now 
closely followed. In fact, our p-meson pole moves gradu- 
ally from the ^ (qq) curve to the (qq) one, as obtained 
in [IS]. However, one must be careful about this conclu- 
sion, since we are disregarding medium related broaden- 
ing which might change the scaling picture observed here. 
Strictly speaking, our results indicate that if the relevant 
density effects amount only to a scaling of f-^, then one 
gets scaling in the pole mass, which is quite consistent 
with the BR idea. 

There is an aditional interesting feature of our results: 
when density is increased further, tttt bound states (first 
sheet poles) appear just below threshold in both chan- 
nels. As we shall see in Sect. 13.21 this result allows for an 
interpretation in terms of a "molecular" classification of 
resonances, which is completely different in both channels. 
In fact, we note that in the a channel the bound state is 
preceded by a doubling of poles in the second sheet. This 
will be further discussed in Sect. 13.21 The appearance of 
pole doubling and bound states has been also analyzed in 
other works |30l|32l|33l. 



3.1 Comparison with nuclear matter experiments. 

As discussed above, apart from Heavy Ion Collisions, there 
are several dedicated experiments on production reactions 
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Fig. 3. Comparison between the (t//o(600) mass (extracted from the lAM pole) at finite temperature and the root of the 
condensate, for the physical pion mass and near the chiral limit. 
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Fig. 4. JAM poles with varying /tt. The numbers attached to the pole positions indicate .fiT(p)(p/ po in brackets) in MeV 
according to ([7|. The circles and open triangles denote second sheet poles, while black triangles refer to first sheet ones. 




0.5 0.6 0.7 0.8 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

./,(p)/y.(o) mp)/mo) 

Fig. 5. Scaling of the mass, compared with fl{p)/fl{0) ~ (qq) (p) / (qq) (0) and with fn{p)/ f-^{0). In the / = J = case, the 
masses displayed correspond to the lowest masses of the second-sheet poles. 
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in nuclei where our analysis can be useful, namely, those 
related to threshoM enhancement in tttt production and 
the modification of dilepton spectra from p decays. 



The threshold enhancement observed in [23ll24[[25] in 
the I = J = channel is qualitatively reproduced in our 
"/,r scaling" finite-density approach. Comparing with our 
finite-T analysis in the previous sections, our main conclu- 
sion is that, due to the broad nature of the a at finite T, 
threshold enhancement is not visible at finite T and it can 
only take place if chiral restoring finite-density effects are 
strong enough, compared for instance with medium effects 
we have neglected which could also increase the a width, 
preventing the pole from approaching the real axis, as in 
the finite-T case, and erasing the threshold effect. This 
does not seem to be the case when many-body interac- 
tions are properly accounted for |26U29j. which gives sup- 
port to this simple approximation. In fact, the observed 
effect is in reasonable quantitative agreement with previ- 
ous theoretical works [27j.,28,30^ and it is not very strong 
for densities p<po (the range available experimentally) in 
line with the experimental results in |24j . 



As for dileptons in nuclear matter, a parameter directly 
measurable in those experiments is the coefficient of the 
linear density term for the mass of the p(770) meson: 



Mjp) 
M{0) 



1 - a- 



Po 



(8) 



The experimental values available so far are a = 0.092± 
0.002 obtained by the E325-KEK collaboration [M] and 
a — 0.02 ± 0.02 measured by the JLab-CLAS experiment 
[TT] . We have performed a linear fit of our pole results in 
the / = J = 1 channel for = 93,85,80,75 MeV, i.e, 
up to p<po to be closest to the experimental situation. 
Our fit gives a — 0.2, which is closest to the result in [14] 
than to |17| . although a bit above the experimental value. 
Our result is also in agreement with Brown-Rho scaling 
[l5] and QCD sum rules [16] which predict a ~ 0.1 — 0.2. 
It is important to remark again that so far we ignore all 
medium broadening effects, which, unlike the case of the 
CT, may be crucial in this emphasized in many- 

body treatments [SUTI ITUIITB] . In fact, in [33] it was reahzed 
that QCD sum rules themselves do not provide a unique 
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constraint for the in-medium mass and width variation, 
unless additional model assumptions are made. In partic- 
ular, it is shown that if one assumes that the width is 
not increased then automatically the mass drops. This is 
the scenario we have recovered in our present approach, 
which, as commented several times above, does not mean 
that this is the physically relevant case. 

3.2 "Molecular" classification of resonances. 

One of the objectives of this work is to analyze the qq 
structure of light meson resonances from a thermal and 
finite-density viewpoint. We have already discussed in Sec- 
tion [2?2] that our results for the thermal a are not consis- 
tent with its pure qq nature. 

The situation changes qualitatively with the T = 
finite-density dependence obtained in Section [3l We see in 
Figure[5]that the a pole follows quite well the same pattern 
of a 0(4) (T-like qq state. One may then wonder about the 
implication of this for the qq nature of the in- medium a. 
We want to point out here that actually one can gain very 
useful information about that by looking at the behavior 
of the poles near threshold. Our argument is supported by 
the classification of resonances lying near threshold given 
in [50j and based on the effective range approximation. 
In those works it is stated that although generally it is 
difficult to extract properties about the "internal" nature 
of resonances from its decay products (scattering poles), 
when the pole lies near threshold, a general rule can be 
applied: a "potential" or "molecular" resonance shows up 
as a single pole near threshold, while for a gq-like state 
two poles near threshold appear very close to another in 
different Riemann sheets. This classification was originally 
applied to states like the /o(980), which lies very close to 
the KK threshold. 

Our claim here is that this classification argument can 
equally be applied to light resonances if medium effects 
drive the poles to the real axis. This is exactly the situa- 
tion when we include density effects only through Eq. ([7]) , 
as it is clear from the results in Figure HI although bear- 
ing in mind that the effective values of for which the 
poles reaches the real axis are too small to trust entirely 
our ChPT-based approach. We see that the a pole follows 
a clear "molecular" pattern, since the pole that remains 
close to threshold and eventually becomes a bound state is 
well separated from the second-sheet pole that lies below 
threshold. We interpret this as a coexistence of two states 
for high enough densities: a tttt "molecule" and a virtual 
state which behaves as a "chiral partner" of the pion, in 
the sense that it tends to become degenerate in mass with 
it (albeit with different quantum numbers) following the 
order parameter. Note that the pole we represent in Figure 
[5] is precisely this virtual state. This picture is in contrast 
with what we observe for the p channel, where two nearby 
poles remain below threshold, one of them moving to the 
first sheet and becoming a bound state. This is clearly 
a qq, non-molecular scenario, according to the previously 
discussed classification. We remark that a similar picture 
for the real axis poles has been obtained in |30j . where 



the density dependence is also parametrized in and in 
vacuum by increasing the quark mass in order to compare 
with lattice results [5Tj . 



4 Finite temperature and density in a 
many-body unitarized approach. 

In this Section we extend our previous results in the / = 
J = channel by incorporating finite nuclear density and 
temperature effects in a many-body description of pion 
dynamics. We follow the line of Refs. 26,31J where the 
effective tttt scattering amplitude in cold nuclear matter 
was evaluated in a chiral unitary framework. For technical 
reasons which we discuss below, it is convenient to use a 
different unitarization scheme for the t-matrix, which how- 
ever provides similar results in vacuum to those from the 
lAM, i.e., the light meson-meson resonances are dynami- 
cally generated in the scalar channel and many scattering 
observables are described in good agreement with experi- 
ment [5^1551151 ■ 

The idea is to solve the coupled-channel Bethe-Salpeter 
(BS) equation for the partial wave scattering amplitude 
(in matrix notation), 

T = V + VGT, (9) 

where the potential (kernel) V of the equation is provided 
by the lowest order tree level amplitude from the chiral 
Lagrangian. In Eq. 

' J {2t:Y {P^ -q^Y -q-^ -ml + ie 

^ {q°f - - mf + is ^^^"^ 

stands for the intermediate two-particle meson-meson Green's 
function of channel i [G is diagonal), where P — (P°, 0) is 
the total external momentum in the center of mass frame 
of the two pions (rest frame with respect to the nuclear 
medium), with s = {P^Y . In principle, both V and T en- 
ter off-shell under the momentum integration {VGT term) 
of the meson-meson loop. However, as it was shown in 
Ref. [52], the (divergent) off-shell contributions of V and T 
in the s-wave interaction can be reabsorbed in a renormal- 
ization of the bare coupling constants and masses order by 
order. Therefore, both V and T can be factorized on-shell, 
leaving the four-momentum integration only in the two- 
particle meson- meson propagator, cf. Eq. (jlOp . An alterna- 
tive justification for solving the Bethe-Salpeter equation 
with on-shell amplitudes can be found in the framework 
of the N/D method, applied for meson- meson [55] and 
meson-baryon [56] interactions. We are thus left with a 
set of linear algebraic equations with trivial solution, 

T= [l-VG]-^V . (11) 

Since we are interested in the a meson energy region (and 
particularly close to the two-pion threshold) , for simplicity 
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^ N ^ ^ * ' 



Fig. 8. Diagrammatic representation of the Bethe-Salpeter ' (a) ' ' (b) 

equation for tttt scattering. r> n- ^ j i j- • i-t • 

^ ° big. 9. Pion tadpole diagrams m tttt scattermg 



we shall work in a single-channel approach (the coupled 
KK channel lies far above in energy and has little effect at 
low energies). The conditions under which this unitariza- 
tion procedure is equivalent to the lAM were discussed 
in [53, 54]. The main difference of the present approach 
as compared to the full I AM amplitude lies in the 0{p'^) 
contribution (which here comes from the s-channel meson- 
meson loop and no tadpole or C4 tree-level terms are in- 
cluded) and the absence of t- and u-channel diagrams (or, 
in other words, the left-hand analytical cut is ignored). 
Still, the present scheme dynamically generates the a pole 
with similar properties as in the lAM and the experimen- 
tal phase shifts in tttt scattering are well reproduced [52l 
[531154] . Also note that Eqs. (Il|4|) are equally satisfied in 
this approach. We show in Fig. [8] a diagrammatic rep- 
resentation of the series of s-channel diagrams which is 
summed in the BS equation. 

The calculation of the thermal amplitude proceeds by 
first reanalyzing the on-shell factorization of V and T 
in the BS equation. At finite temperature, the divergent 
contribution from off-shell terms in the one-loop ampli- 
tude has the same structure as in vacuum, but a finite, 
temperature-dependent part survives which cannot be cast 
as a (vacuum) renormalization of /^r or ttIt^. This contri- 
bution can be accounted for as a temperature correction 
to the O(p^) kernel. 



3/^ 



Vonis) + 



3fS 



[12) 



with 7^ = (277^) ^ dqq'^nB{(^q)/ujq, Von{s) = ~{s - 
mJ/2)//^ and cUg = + (we follow in this section 
the normalization of partial waves given in [521l53j , which 
differs from the one in Sect. [2] in a factor IBtt). In order to 
keep as close as possible to the physics described by the 
I AM amplitude, we have also considered the finite-T con- 
tribution from tadpole terms, which becomes relevant as 
the a meson pole is driven towards the two-pion threshold 
and its behavior is no longer dominated by the two-pion 
phase space j42] . The corresponding diagrams are depicted 
in Fig. [9] and, to the lowest order in the chiral counting, 
they emerge from 0{p^) interaction terms with up to six 
meson fields. The finite-T correction to the tree level am- 
plitude to account for these terms reads 



'5Kld(s) 



with I2 = m 



25 , 



(13) 



2 tT 

TT ^0 



In Fig. [TO] (left) we show the a pole trajectory for the 
thermal calculation in the BS equation approach with the 
0{p^) kernel including the thermal corrections discussed 



above. As we can see, the evolution of the pole follows 
quite a similar pattern as in the JAM, although the posi- 
tion of the (7 mass reaches values further below the two- 
pion threshold for high temperatures. As in the lAM, the 
a pole stays far from the real axis, indicating a substan- 
tial width for temperatures up to about 200 MeV, despite 
the low value of Mp. The stronger attraction felt by the 
a meson in this approach as compared to the lAM seems 
to reflect that missing thermal contributions from t- and 
M-channel dynamics are relevant in this energy region and 
provide a repulsive contribution. Nevertheless, we obtain 
a qualitatively similar physical behavior of the a pole in 
both the BS equation approach and the I AM. Thus, we 
shall use the former as a starting point to incorporate 
finite density effects. As a further test, we have also calcu- 
lated the a pole trajectory in the effective "/^ scaling" sce- 
nario, which we depict in the right hand panel of Fig. 1101 
We obtain a similar result as in the lAM regarding the a 
pole collapsing onto the real energy axis for values be- 
low approximately 50 MeV. We also find a pole-doubling 
effect which follows the "molecular" pattern discussed in 
Sect. [3] Differences are observed at the numerical level as 
the close-to-threshold behavior of the a pole is dictated 
by the relative weight of the 0{p^) and 0{p'^) amplitudes, 
which are different for the two approaches discussed in this 
work. Finally, we have also studied threshold enhancement 
in the tttt amplitude for decreasing values of /^r , which we 
omit here as our results resemble very much those depicted 
in Figs. [5] and [7] for the lAM calculation. 

The introduction of nuclear density effects on top of 
the temperature follows by a renormalization of the pion 
propagator in the hot and dense medium. In cold nu- 
clear matter, the pion spectral function exhibits a mixture 
of the pion quasi-particle mode and p-wave particle-hole 
(ph), Delta-hole (Ah) excitations f57^. The lowest order, 
irreducible p-wave pion self-energy due to ph and Ah ex- 
citations reads 

^^nn-^+ttan-^ (90, g; T) = 
—] [UNN-^{qo,q;T) + u^N-^iqo,q;T)] , 

nirr J 

(14) 

where U stands for the finite temperature Lindhard func- 
tion, which we evaluate in Imaginary Time Formalism 
(ITF) [58l[38] , and the density dependence enters U through 
the baryon chemical potential. We use phenomenological 
ttNN and ■nN A coupling constants determined from anal- 
ysis of pion nucleon and pion nucleus reactions, ji<s Im^^ = 
0.007244 MeV-i and /zi/Zw = 2.13. The strength of the 
considered collective modes is modified by repulsive, spin- 
isospin and N A short range correlations ^57], which 
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Fig. 10. Left: Temperature dependence of the a complex pole in the BS equation approach including thermal corrections to 
the 0{p^) kernel (the temperature interval between each point is 5 MeV). Right: Pole trajectory from scaling. The numbers 
attached to each point indicate the value of /,r in MeV. 



we include in a phenomenological way with a single Landau- 
Migdal interaction parameter, g' = 0.7. The RPA-summed 
(retarded) pion self-energy then reads 

7TP((?o,<7;T) = 

(^) ' ^-(g ') 9 ' [Unn-^ {qo,q; T) + C/^at-i (go, q; T)] 



1 - 



(It) 9' [UnnMio^Q^T) + U^M-^iqo,q;T)] 



(15) 



where we have accounted for the finite size of ttNN and 
ttNA vertices with the hadronic monopole form factors 
F^iq"^) = Al/{Al -q^), with = 1300 MeV. 

The in-medium pion propagator modifies the analyt- 
ical structure of the meson-meson loop function, G. At 
lowest order in a density expansion (number of baryon- 
hole irreducible insertions) the in-medium irn amplitude 
at one loop includes diagram (a) in Fig. 111! in which one 
of the intermediate pions excites a ph (or Ah) bubble (and 
similarly with the lower pion line), on top of the vacuum 
TTTT loop (Fig. [m second diagram). Again one is obHged 
to check whether the on-shell factorization of V and T is 
still valid in a nuclear medium. This is actually the case 
for TTTT scattering in the scalar channel as it was shown in 
[25] . The argument is as follows: in addition to diagram 
(a), chiral symmetry requires a set of meson baryon con- 
tact terms, depicted as diagrams (b-d) (which can be seen 
as contributions to the wN — > ttttA^ amplitude by cutting 
simultaneously the ph bubble and the lower pion line in 
diagrams (a-d)). It turns out that the contribution from 
the off-shell part of V in diagram (a) exactly cancels with 
the sum of the amplitudes from diagrams (b-d), leaving 
us with diagram (a) with each of the tttt vertices factor- 
ized on-shell. Therefore, the algebraic solution of Eq. ([9]) 
is no altered and one only has to replace the vacuum pion 
propagators in G, cf. Eq. (jlOp . by the in-medium ones. 



X X + )< + + ^ — y + 



Fig. 11. Lowest order baryon-hole insertions in the meson 
meson loop function including vertex corrections. 



temperature as we have checked explicitly, in ITF, that 
the same cancellation of off-shell terms takes place in the 
thermal amplitude. 

In order to re-evaluate G in the hot and dense medium 
we use the spectral (Lehmann) representation of the pion 
propagator. The final expression for G, once continued 
onto the real energy axis, reads 



GiP;p,T) 
2^ 



1 



1 



with 



po - f2 + ie + [2 + ie 



Fin) 



(16) 



,0\2 



TTtt] . The argument holds at finite 



X S^iE+,q)S^iE^,q) , (17) 

where S"^ = —tt^^IirDt^ is the spectral function of the 
retarded pion propagator and E± = (i? ± u)/2 (we have 
omitted here the contribution from diffusion poles which 
typically provides a small correction in the time-like re- 
gion). Note that Im G(P") ^-F{P")/2 and thus F plays 
the role of a generalized in-medium two-pion phase space 
including both temperature and nuclear density effects. In 
the pure thermal case, F{y/s) = axis) /Sir 9{s — 4m^). 

In Fig. [T2]we show the phase space function, i^(P"), 
at normal nuclear matter density {po — 0.17 fm^"^) and 
different temperatures in the range T = — 150 MeV. We 
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Fig. 12. Two-pion phase-space function at finite nuclear den- 
sity and temperature. 



0( p) 



Fig. 13. 0{l,p,p^) contributions to the meson meson loop 
function from irreducible pion selfenergy insertions. 



have also depicted the vacuum and thermal cases for com- 
parison. The very first difference that one observes at finite 
nuclear density is the appearance of strength below the 
two-pion threshold, which is absent in the thermal case. 
It can be understood on the basis of the baryon-related 
interaction mechanisms of the pion discussed above. At 
small nuclear densities, one of the intermediate pions may 
excite a ph pair whereas the other is placed on the mass 
shell (see the second diagram in Fig. [T5)) . This mechanism 
is responsible for the strength right below 2m7r and low- 
ers the threshold down to mjr. For increasing density, the 
probability for the two pions to be absorbed by baryon- 
hole excitations sets in (third diagram in Fig. [T3|) . which 
builds additional strength below and shifts the thresh- 
old practically down to P° = (note the smaller size of 
this 0{p^) contribution with respect to the excitation of 
one single ph bubble). Beyond P° = 2m^, F(P°) also 
exhibits a remarkable enhancement with respect to the 
vacuum and thermal cases, indicating an increased phase 
space for a — > tttt decays. This reflects the widely spread 
structure of the pion spectral function in the medium and 
the considerable attraction experienced by the pion quasi- 
particle peak. 

When the temperature is increased, the low energy re- 
gion looses some strength, as the phase space for NN~^ 
excitations is smeared off by the thermal motion of the 
nucleons. Above = ttt-tt, Bose enhancement is more ef- 
fective as one of the pions is placed on-shell and one can 
appreciate some increase over the T = result. The lat- 
ter effect is strongly magnified right above the two-pion 



threshold, since (i) the pion spectral function is popu- 
lated at low energies and (ii) the pion quasiparticle peak 
is strongly attracted in the medium. Note that the rapid 
energy dependence of the p-wave pion selfenergy is also 
responsible from the quick increase of phase space right 
beyond the opening of the tttt channel. 

Wc finally show in Fig.[T3]the imaginary part of the tttt 
amplitude for different densities and temperatures, which 
we can compare to the thermal calculation and the 
scaling" scenario in order to study threshold enhancement 
effects. We recall here that in the vacuum and thermal 
cases ImT vanishes below the two-pion threshold. Our re- 
sult, from a dynamical calculation of many-body effects 
in the intermediate two-pion state, exhibits a sizable ac- 
cumulation of strength at and below P*^ = 2m7r already 
at nuclear matter density and zero temperature. On the 
other hand, the amplitude is strongly depleted beyond 
a/s — 400 MeV as compared to the vacuum case, an ef- 
fect that is not shared by the reduced- calculation at 
the equivalent density (cf. Fig. [7|) which does not account 
for baryon-related cr-decay mechanisms. Only for much 
smaller values of the strength focuses at the vacuum 
threshold, as the a pole lies very close to the real en- 
ergy axis; the resulting enhancement in the amplitude in 
that case is by far larger than obtained here, which is re- 
lated to, but not fully driven by the a pole behavior as we 
explain below. With respect to the thermal calculation, 
finite density additionally softens the amplitude at high 
energies and brings strength below the vacuum threshold, 
a feature which is linked to the presence of a baryonic 
medium and the excitation of pionic collective modes, as 
discussed above. 

We have also studied the behavior of the a pole at 
finite nuclear density. The analytic continuation of T in 
this case is not trivial and we have used an approximated 
prescription, namely to add to the evaluation of G in 
the 1st Riemann sheet the discontinuity on the real axis 
with respect to the unphysical (2nd) sheet, DiscG(P°) = 
— 2iImG(P°). Still, this procedure provides the qualita- 
tive in-medium behavior of the a pole [53], although its 
absolute position in the complex plane carries some uncer- 
tainty (which we believe is superseded by other theoretical 
uncertainties in the evaluation of the pion selfenergy) . The 
combined effects of temperature and nuclear density accel- 
erate the migration of the cr pole towards the tttt threshold. 
For instance, a.t p = po/2, the a mass reaches Mp = 2m^ 
at about T ~150 MeV, whereas at normal nuclear density 
Mp ~ 300 MeV already at T = and it quickly reaches the 
two-pion threshold at about T — 100 MeV. This is corre- 
lated with the cusp structures observed in the amplitude 
at threshold, cf. Fig. [TH however, the strength observed 
at lower energies is linked to the many-body pion dynam- 
ics. In spite of this, the a remains as a broad resonance at 
nuclear matter density and temperatures approaching the 
transition one, similarly to what happens in the thermal 
calculation (at comparatively higher temperatures) and at 
variance with the simplified "/^ scaling" approximation. 
In fact, if we keep increasing the density, at some point 
the a pole crosses below 2m7r (one should not trust our 
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implementation of medium effects far beyond p — po, but 
as an exercise it provides information about the phase- 
space behavior of the resonance). For instance, at p = 2po 
and T = 100 MeV the a mass from the pole lies about 
70 MeV below the two-pion threshold but still we find 
Fp ~ 150 MeV. A more detailed investigation of possible 
pole-doubling effects in this approach is on-going and will 
be reported in a future work. 



5 Conclusions. 

We have presented an analysis of the behavior of tttt scat- 
tering amplitudes in Unitarized Chiral Perturbation The- 
ory with medium effects incorporated in several ways. In 
particular, we have been focused on the behavior with fi- 
nite temperature and nuclear density of the p(770) and 
/o(600)/(j resonances, which are generated dynamically 
within the Inverse Amplitude Method (lAM). 

By considering only thermal effects on the lAM tttt 
amplitudes, the p exhibits a considerable broadening with 
a small mass decrease as temperature increases, whereas 
the a mainly decreases its mass, effectively signaling chi- 
ral symmetry restoration, although it still remains as a 
broad resonance even at the transition temperature. The 
broadening obtained in our approach for the p meson at fi- 
nite temperature is compatible with the spectral function 
analysis from dilepton spectra in the recent experiment 
by the NA60 Collaboration. The evolution of the p mass 
with temperature does not scale as the condensate, which 
renders our results in quantitative disagreement with the 
Brown-Rho scaling scenario. The fact that the a pole re- 
mains far from the real axis even at the two-pion threshold 
when only temperature effects are considered implies no 
significant threshold enhancement for the scattering am- 
plitude, which has been advocated as a precursor of chiral 
symmetry restoration. We neither observe a scaling of the 
a mass with the quark condensate, which indicates that 
the /o(600) / <T resonance dynamically generated in our uni- 



tarized chiral approach has a non-qq component which is 
relevant near the phase transition. 

By introducing finite nuclear density the picture changes 
dramatically. In a first approximation we have incorpo- 
rated the effect of a nuclear medium by decreasing /^r ac- 
cording to the GOR relation to linear order in density. At 
sufficiently low (high) values of /^r (density), the p and 
<j poles collapse onto the real energy axis at the thresh- 
old energy, which is preceded by a significant threshold 
enhancement in the scattering amplitudes. We have dis- 
cussed these effects in the context of recent results from 
resonance production in finite nuclei and our results are in 
line with the experimental observations. A detailed analy- 
sis reveals that when the resonance pole is close to the real 
axis it splits into two states in separated Riemann sheets. 
The markedly different properties of these double poles 
for the p and a channels allows us to classify these tttt 
resonances according to their internal structure: whereas 
the p meson presents a clear predominant qq behavior 
for high densities (the two poles stay close to threshold, 
one of them migrating to the 1st Riemann sheet as a tttt 
bound state), the a exhibits a "molecular" behavior (one 
of the poles stays close to threshold, well separated from 
the other one which evolves to lower energies to become 
degenerated with the pion). The mass scaling from the cr 
and p pole with /^r {p) follows the quark condensate evolu- 
tion and therefore is compatible with a Brown-Rho scaling 
scenario, although one should keep in mind that relevant 
finite density mechanisms are neglected in this approxi- 
mation. 

Finally, we have improved our implementation of finite 
nuclear density (and temperature) effects by considering 
a microscopic calculation of many-body pion dynamics in 
TTTT scattering. We have chosen a different unitarization 
scheme for the tttt scattering amplitude, namely to solve 
the Bethe-Salpeter equation for the lowest order ChPT in- 
teraction. Despite differences in the amplitudes at ©(p^), 
this scheme essentially provides the same results as the 
1AM and allows a systematic analysis and resummation of 
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a relevant class of pion interaction mechanisms with the 
nuclear medium. The pion interactions with the medium 
are encoded in the single-particle pion sclfencrgy, which 
accounts for the excitation of p-wave particle-hole and 
Delta-hole components as well as short distance correla- 
tion effects from nucleon-nucleon and Delta-nucleon in- 
teractions. The opening of baryon-related channels on top 
of cr ^ TTTT at finite density extends the available phase 
space to lower energies and therefore the tttt scattering 
amplitude exhibits an increased strength at and below the 
two-pion threshold, which is magnified at finite tempera- 
ture as a consequence of Bose enhancement on the low 
energy modes of the tttt intermediate states. Such an ef- 
fect has been found to provide a satisfactory description 
of the data from the two-pion photoproduction reaction in 
nuclei when comparing the mass spectrum in the neutral- 
vs charged-pion channels for different nuclei [29ll25] , where 
nuclear densities of the order of po and below are explored. 
In our analysis we have considered both finite temperature 
and nuclear density, thus extending the applicability of the 
present approach to other experimental scenarios such as 
the forthcoming Heavy-Ion physics program at FAIR. As 
compared to the purely thermal calculation, the attrac- 
tive interaction mechanisms of the pion at finite density 
accelerate the migration of the a pole towards the two- 
pion threshold. The threshold enhancement observed in 
the TTTT amplitude is correlated to the evolution of the a 
pole towards (and below) 2m,r- However, differently from 
the reduced- /,r result, the tr pole stays far from the real 
axis indicating a sizable decay width of the resonance at 
densities as high as 2pQ and temperatures close to the 
transition one. 

As a continuation of this work we plan to implement a 
similar many-body analysis of finite density effects for the 
electromagnetic pion vector form factor, an extension of 
the present work to the SU{Z) case where other resonances 
and heavier meson states come into play, and to introduce 
the effect of a finite meson-number chemical potential. We 
will report on these studies elsewhere. 
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